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Abstract 

We show that the minimal speed for the existence of monotonic fronts of the equation 
Ut = (u m ) xx + f(u) with /(0) = f(l) = 0, m > 1 and / > in (0, 1), derives from a 
variational principle. The variational principle allows to calculate, in principle, the exact 
speed for general /. The case m = 1 when f'(0) = is included as an extension of the 
results. 
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Several problems arising in population growth [3], ||, combustion theory || U, chemical 
kinetics ||, and others ||, lead to an equation of the form 

where the source term F(p) represents net growth and saturation processes. The flux j is 
given by Fick's law 

where the diffusion coefficient D(p) may depend on the density or in simple cases be taken 
as a constant. In one dimension this leads to the equation 

In what follows we shall assume that 

F(p)>0 in (0,1), and F(0) = F(l) = 0, (16) 

restrictions which are satisfied by several models. When the diffusion coefficient is constant 
and the additional requirement F'(0) > is satisfied, the asymptotic speed of propagation 
of localized small perturbations to the unstable state u = is bounded below and in some 



cases coincides [[?[] with the value cl = 2 J F'(0) which is obtained from considerations on the 
linearized equation ||. However, when either F'(0) = or D(p) is not a constant, no hint for 
the speed of propagation of disturbances can be obtained from linear theory alone. A common 
choice for the diffusion coefficient is a power law, case with which we shall be concerned here. 
Therefore the equation that we study is 

% = (P m ) xx + F(p) (2a) 

with 

F(0) = F(l) = 0, and F > in (0, 1). (2b) 

Aronson and Weinberger j7|, ^| have shown that the asymptotic speed of propagation of 
disturbances from rest is the minimal speed c*(m) for which there exist monotonic travelling 
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fronts p(x,t) = q(x — ct) joining q — 1 to q — 0. The equation satisfied by the travelling 
fronts is 

(q m )z Z + cq z + F(q) = (3a) 

with 

q(— oo) = 1, q > 0, q < in (— oo,co>), g(co>) = (36) 

where z = x — ct. The wave of minimal speed is sharp, that is, ui < oo when m> 1 §. 

An explicit solution is known Q for the case F(q) = q(l — q) and m = 2, the waveform 
is given by 

q{z) 



1 - 
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J + 

and it travels with speed c*(2) = 1 (here [x] + = max(x, 0)). Recently the derivative dc/dm at 
m = 2 has been calculated by two different methods. Its value is -7/24 [§, |i~0| . Other exact 
solutions for different choices for m and F have been given in ||11|| . 

The purpose of this work is to give a variational characterization of the minimal speed 
c*(m) for Eq.(3) when m > 1, and as a byproduct for the case m = 1 when F'(0) = 0, both, 
cases for which no information is obtained from linear theory. The case m = 1 with -F'(O) > 
has been studied elsewhere ||13|| . Lower bounds have been obtained on the minimal speed 
12| ; the present results allow its exact calculation for arbitrary _/. 
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Since the selected speed corresponds to that of a decreasing monotonic front, we may 
consider the dependence of its derivative dq/dz on q. Calling p(q) = —q m ~ 1 dq/dz, where the 
minus sign is included so that p is positive, we find that the monotonic fronts are solutions of 

P^--p+-q m ~ 1 F(q) = (4a) 
dq m m 

with 

p(0)=p(l)=0, p>0 in (0,1). (46) 

Although the wave of minimal speed is sharp and therefore g'(0) < 0, by its definition p(0) = 
is true. We now show that the minimal speed c*(m) follows from a variational principle whose 
Euler equation is Eq.(4a). 
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Let g be a positive function such that h = —g' > 0. Multiplying Eq.(4a) by g/p and 
integrating we obtain after integration by parts, 

/ 1 [^ m - 1 F(g)^ + %Mg)] dq 



c 



By Schwarz's inequality, since, q, F, g and h are positive we know 



(5) 



]_ q m-i F 9_ + hp > 2J—q m - 1 Fgh (6) 
m p V m 



and therefore, replacing in Eq.(5) we have 



c> 2 IoV m Q m l Fghdq _ 
Jo 1 9dq 



This bound has been already given by us fll2| . We now show that it is always possible to 
find a g(q) such that the equality in Eq.(6) and therefore also in Eq.(7) holds. We do so by 
explicit construction of such a function g. The equality in Eq.(6) holds if 



]_ q m~l F 9 = hp ^ 

m p 



Let v(q) be the positive solution of 



and choose 



v' c 



(9a) 
rap 



9 = - (96) 
v' 



We have then 



v" (v') 2 c 



2 P' = —q m - l F(q) 
v v mp z m p^ 



where we have used Eq.(9a) to eliminate v' and Eq.(4a) to eliminate p' . Therefore, 

h = -g' = /- = — 2 q m - l Fg > (9c) 
[y'y mp z 

where we have made use of Eqs.(9a) and (9b). With this expression for h, we can see that 
Eq.(8) is satisfied. In addition we must check that g as we have defined it is such that its 
integral exists. In fact as it exists and moreover one can always normalize g so that g(0) = 1 
and g(l) = 0. From the definition of g we obtain 



mp(q) 
g(q) = exp 



go rap 



where < go < 1- Since p(l) = and p is positive between and 1 it follows that g(l) = 0. 
At zero no divergence occurs, as we now show. Call c = c/m and f(q) = q m ~ l F(q)/m. Then 
Eq.(4a) reads 

pp' — cp + f = (10a) 

with 

/(0) = /(l)=0 and /'(0) = 0. (106) 

For this case Aronson and Weinberger have shown that p(q) approaches the fixed point 
q = as p = cq = cq/m. Then, near 0, v'/v ~ \ jq or v ~ q and from its definition g(0) = 1. 
Then the integral of g exists. We have shown then 

*, x n IoVrnq m -^Fghdq 

c {mj = max 2 — j . (11) 

Jo gdq 

where the maximum is taken over all functions g such that 

0(0) = 1, 0(1) = and h = -g' > 0. 

It is perhaps of some interest to verify explicitly that the Euler equation for the maximizing 
g is indeed Eq.(4a). Let us study the maximization of the functional 

J m (g) =2 J s/mq^Fghdq 

where h = —g' > subject to 

/ g(q)dq=l. 
Jo 

The Euler equation for this problem is 



\mq m - x Fh d / mq m - 1 Fg\ 

A + V — i — + ^(,v — ^ — j =0 

where A is the Lagrange multiplier. Using the expression given in Eq.(9c) for h we see that 
this is exactly Eq.(4a) with the Lagrange multiplier A = — c. 

As an application we shall consider the case F(q) = q(l — q) and m = 2 for which the 
exact solution is known. Take as the trial function g(q) = (1 — q) 2 . Then we obtain 

C>4 AV'-^ = 1 

" /o(l-?) 2 <<? 



the exact value, which shows that this is the function g for which the maximum is attained. In 
addition, due to the existence of the variational principle we may use the Feynman-Hellman 
formula to calculate the dependence of c(m) on parameters of F. We illustrate this by applying 
it to the calculation of dc/dm at m = 2. Taking the derivative of Eq.(lO) with respect to m 
we obtain 

Evaluating at m = 2, with g(q) = (1 — q) 2 we obtain 

*(2) = 3 J f 9 (l- 9 )'(l + 21og,)*=4 

the value previously obtained by other methods. 

A fast estimation of the speed for other values of m can be obtained with simple trial 
functions. In Fig. 1 we show lower bounds for F = q(l—q) using as trial functions g\ = (1— q) 2 
and #2 — (1 — With the first trial function we have the exact value at m = 2. The dotted 
line is the line of slope -7/24 that coincides with the tangent at m = 2. For larger m a better 
estimate is obtained using g^. The dashed line is the curve \Jl/m which has been suggested 
by Newman |l[ as the best fit to his numerical results. With better choice of trial functions 
the exact value can be approached arbitrarily close. 

Finally we observe that the case m = 1 when F'(0) = follows directly here. Repeating 

the procedure starting now from equation (10), one obtains, 

o loVFghdq 

c = max 2 t 

Jo 9dq 

where the maximum is taken over all functions g such that 

#(0) = 1, 0(1) = and h = -g' > 0. 

To show this we have used v'/v = c/p and g = 1/v' and the asymptotic behavior described 
above. 
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